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Abstract 

Climate models predict a range of changes in tropical forest regions, including increased average temperatures, decreased 
total precipitation, reduced soil moisture and alterations in seasonal climate variations. These changes are directly related to 
the increase in anthropogenic greenhouse gas concentrations, primarily C0 2 . Assessing seasonal forest growth responses to 
climate is of utmost importance because woody tissues, produced by photosynthesis from atmospheric C0 2 , water and 
light, constitute the main component of carbon sequestration in the forest ecosystem. In this paper, we combine intra- 
annual tree growth measurements from published tree growth data and the corresponding monthly climate data for 25 
pan-tropical forest sites. This meta-analysis is designed to find the shared climate drivers of tree growth and their relative 
importance across pan-tropical forests in order to improve carbon uptake models in a global change context. Tree growth 
reveals significant intra-annual seasonality at seasonally dry sites or in wet tropical forests. Of the overall variation in tree 
growth, 28.7% was explained by the site effect, i.e. the tree growth average per site. The best predictive model included four 
climate variables: precipitation, solar radiation (estimated with extrasolar radiation reaching the atmosphere), temperature 
amplitude and relative soil water content. This model explained more than 50% of the tree growth variations across tropical 
forests. Precipitation and solar radiation are the main seasonal drivers of tree growth, causing 19.8% and 16.3% of the tree 
growth variations. Both have a significant positive association with tree growth. These findings suggest that forest 
productivity due to tropical tree growth will be reduced in the future if climate extremes, such as droughts, become more 
frequent. 
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Introduction 

Tropical forests are being threatened on an unprecedented scale 
by global changes. Temperatures across tropical forest regions are 
currently increasing [1] and are expected to continue to increase 
with a concomitant decrease in precipitation over the next decades 
[2-4] . Climate models predict a range of changes in tropical forest 
regions, including increased frequency of extreme climatic events, 
increased average temperatures, increased atmospheric C0 2 and 
changes in seasonal distribution and interannual variability of 
rainfall [5-9]. Tropical forests play an important role in the 
mitigation of anthropogenic atmospheric C0 2 emissions by 
constituting a major reservoir of terrestrial carbon and a large 
and persistent carbon sink [10-12]. Feedback between tropical 
forests and the local and regional climate has also been 
demonstrated [4,13]. 



Tree growth is linked with atmospheric C0 2 through photo- 
synthesis. The last 20 years have seen a substantial increase in the 
number of publications focusing on the effects of climate on 
tropical tree growth. A search performed on Web of Science in 
March 2013 using the keywords 'climate', 'tropical forest', 
'growth' and 'trees' returned fewer than 15 articles per year 
before 2000 and more than 60 articles in 2012, for a total of 541 
articles focusing on the effect of climate on tropical tree growth. 
Long-term variations in tree growth have been reported on long- 
term forest plots, but the determinants of these variations are still 
being discussed [14-19]. Due to the annual or multi-annual census 
frequency of long-term forest plots, most studies focus on the 
annual or multi-annual variation in tree growth even though most 
tropical forests undergo an intrannual seasonality in climate 
[1,6,20-22]. In single-site-based studies, seasonal rhythms of tree 
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Table 1. Expected tropical tree growth response to climate variables. 





variable 


predicted effect' 


references 


process* 


REW 


+ 


[28,62] 


photosynthesis, xylem tension, stomatal closure, leaf flush 




rainfall 


+ 


[19, 23, 24, 35, 58, 59, 61, 97-90] 


photosynthesis, xylem tension, stomatal closure, leaf flush 






[89,91] 




T mean 




[31-33,92,93] 


photosynthesis kinetic, stomatal closure 


T min 




[14,17,35,66] 


photosynthesis kinetic, stomatal closure 




no 


[94,95] 




T max 




[19,31,33] 


photosynthesis kinetic, stomatal closure 


+ [96] 




no 


[94,95] 




VPD 


no 


[97,98] 


stomatal closure, transpiration 


irradiance 


+ 


[17,21,29,30,64,70,74] 


photosynthesis, phenology 






[21] 






no 


[17,63,70] 




U* 


+ 


[99] 


photosynthesis, transpiration 


a : expected growth response to the climate variable: (+) trees are expected to grow faster with high values of the climate variable, ( — ) trees are expected to grow slower 
with high values of the climate variable. b : biological processes involve in the tree growth response to a given climate variation. VPD is vapour pressure deficit, and 
Friction velocity (L/*) is a climate variable provided by eddy flux data, which is correlated with wind speed. Relative extractable water (REW), is a daily value between 0 



and 1; when REW = 1, the amount of extractable water by the tree is at its maximum and when REW = 0, no water is available for the trees [28]. 
doi:l 0.1 371 /joumal.pone.0092337.t001 



growth have been linked to seasonal variations in water 
availability, rainfall, temperature and solar radiation (Table 1). 

Rain or lack of rain is often implicitly viewed as the main drivers 
of forest dynamics [23], as annual net primary production (NPP) 
positively correlates with the annual sum of precipitation at large 
scales [24] and rainfall seasonality plays a key role in the forest's 
response to climate variability [25]. The relation between the 
amount of rainfall and water availability for trees is not 
straightforward and is determined by various soil and plant 
characteristics (i.e. permanent wilting point, field capacity, root 
distribution). Consequently, water stresses are increasingly esti- 
mated using soil water balance models [26,27], including some 
that are explicitiy designed for tropical forests [28]. Irradiance is 
directly linked to plant photosynthetic capacity, which in turn 



drives carbon uptake and plant growth [29]. The occurrence of 
dry periods linked to cloud-cover reduction was found to enhance 
canopy photosynthetic capacity by 25% throughout Amazonia 
[30]. The effects of rising temperatures on the physiology of 
tropical trees are currently debated within the scientific commu- 
nity [21,32]. Some studies suggest that reductions in photosyn- 
thetic rates at temperatures above 30°C are driven by reductions 
in stomatal conductance in response to higher leaf-to-air vapour 
pressure deficits [31] or by a direct down-regulation of biochem- 
ical processes during CO z fixation [33,33]. Recent studies, 
however, suggest that tropical tree mortality may increase 
significandy with increasing night-time temperatures, while tree 
growth appears surprisingly sensitive to variations in mean annual 
night-time temperatures of 1— 2°C [35]. 



Table 2. Climate datasets used to model seasonal tropical tree growth at a pantropical scale. 



climate variable 


full name 


unit 


origin 


pre 


precipitation 


mm 


CRU-TS3.10.01 


soF 


Extraterrestrial solar radiation 


mm/month as equivalent of evaporation 


CGIAR-CSI 


dd 


cloud cover 


% 


CRU-TS3.1 


Pre 


Potential Evapotranspiration 


mm 


CRU-TS3.1 


swe 


relative soil water content 


unit less 


CGIAR-CSI 


tmp 


mean temperature 


°C 


CRU-TS3.1 


tmn 


minimal temperature 


°C 


CRU-TS3.1 


tmx 


maximal temperature 


°C 


CRU-TS3.1 


vap 


vapour pressure 


hPa 


CRU-TS3.1 


dtr 


temperature amplitude 


°C 


CRU-TS3.1 


a : sol per day (mm/day equivalent) is calculated using the methodology presented in [39] specifically for the 15 th day of each month to describe averages per month. 



Total sol per month (mm/month equivalent) is calculated by multiplying the value of sol for the 15th day of the month by the number of days in the month, 
1 mm.day -1 equivalent of evaporation = 2.45 MJ.m _2 .day _1 . 
doi:1 0.1 371 /joumal.pone.0092337.t002 
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A pantropical analysis of the effects of climate seasonality on 
tropical tree growth is still missing in the literature. Most of our 
knowledge comes from single-site-based studies that often suffer 
from collinearity problems between climate drivers. In this paper, 
we ran a meta-analysis of monthly tropical tree growth at the pan- 
tropical scale in which multiple gradients of climate variables 
allowed us to disentangle the effect of each climate driver on 
tropical tree growth. We focused on the seasonal effect of climate 
on tree growth in tropical forests by using data from papers 
reporting tree growth measurements with a high periodicity (from 
daily to monthly censuses) and global climate datasets. We 
gathered 30 datasets in which growth measurements have been 
recorded for 3412 individual trees from 25 pantropical forest sites. 
This paper has three specific objectives: (i) to find the climate 
drivers of tree growth across tropical forests; (ii) to quantify tree 
growth variations in response to climate among tropical forests; 
and (iii) to give a modeling framework to improve the model 
predictions of seasonal carbon uptake by tropical tree growth in a 
global change context. 

Materials and Methods 

Climate datasets 

We used climate datasets from two sources (Table 2): the 
Climate Research Unit (CRU) at the University of East Anglia 
[36] and the Consortium for Spatial Information website 
(CGIAR-CSI, http://www.cgiar-csi.org). From the CRU, we used 
variables from the CRU-TS3.1 and CRU-TS3. 10.01 monthly 
climate global datasets available at 0.5° resolution from 1901— 
2009: cloud cover {eld, unit:%); precipitation {pre, mm); daily mean, 
minimal and maximal temperatures (respectively tmp, tmn and tmx, 
°C); temperature amplitude {dtr, °C); vapour pressure (vap, hPa); 
and potential evapo transpiration {pet, mm), pet was calculated using 
the grass reference evapotranspiration equation [37,38], which is a 
variant of the Penman-Monteith method using the gridded tmp, 
tmn, tmx, vap and eld. pre was square-root transformed prior to data 
modeling to address heteroscedasticity. From the CGIAR-CSI, we 
used the monthly average of extraterrestrial solar radiation as well 
as the Global Soil- Water Balance [27]. The first dataset defines the 
solar radiation reaching the top of the Earth's atmosphere (sol) and 
is a function of Earth-sun geometry and time of year, sol per day 
(mm/day equivalent) is calculated using the methodology 
presented in [39] specifically for the 15"' day of each month to 
describe averages per month. Total sol per month (mm/month 
equivalent) is calculated by multiplying the value of sol for the 1 5th 
day of the month by the number of days in the month, 
1 mm. day 1 equivalent of evaporation = 2.45 MJ.m 2 .day 
The latter dataset provides hydrological raster data (ESRI Grid 
format) describing actual evapo-transpiration and soil water 
content for a monthly time-series from period 1901-2009 using 
CRU-TS3.1 as the primary climate data input. The monthly time- 
series of relative soil water content swe (unitless, 0-1) is calculated 
as the ratio of the soil water content from the soil water balance 
over the maximum available soil moisture (in mm water per 1 m 
soil depth, from the Digital Soil Map of the World and Derived 
Soil Properties [40,41]) along the ecosystem rooting depth [42]. 

Tree growth data 

A search performed on Web of Science in March 2013 using the 
keywords 'climate', 'tropical forest', 'growth' and 'trees' returned 
541 articles focusing on effect of climate on tropical tree growth. 
Among these publications, our analysis used only the publications 
with seasonal tree growth data (from daily to monthly censuses) 
available in the article graphics, available in an online repository 



4 



March 2014 | Volume 9 | Issue 3 | e92337 



Climate Effects on Seasonal Tropical Tree Growth 




Figure 1. Locations of the 25 study sites and their countries (grey areas). 1: Attapadi; 2: Budongo; 3: CPM; 4: El Palmar; 5: FLONA SFP; 6: 
Guanacaste; 7: Ibicatu; 8: Kakamega; 9: La Barcinera; 10: La Selva; 11: Lamto; 12: Luki forest; 13: Marajoara; 14: Muara Bungo; 15: Munessa-Shashamene 
Forest; 16: Paracou; 17: Pinkwae; 18: RBSF; 19: RFC; 20: Rio Cachoiera; 21: Selangor plantation; 22: SERS; 23: Tapajos; 24: Tinte Bepo; 25: ZF-2. 
doi:10.1371/journal.pone.0092337.g001 



or provided by the authors. The tree growth data were extracted 
from the article graphics using WebPlotDigitizer 2.5 (http:// 
arohatgi.info/WebPlotDigitizer/, Table 3). Three types of data 
were found: (i) single observations by tree and time step at a given 
site (see [43]); (ii) mean and standard deviation (SD) by species at a 
given site for each time step (see [44]); and (iii) only the mean and 
SD of all trees at a site for each time step [35] . In addition, these 
data may come from measurements with dendrometer bands, 
electronic point surveys, tree ring analyses or graduated tapes. The 
minimum diameter at breast height (DBH) used in the articles is 
10 cm. Among sites, Selangor and Muara Bungo are forest 
plantations. Finally, they have different time step between 
measurements: monthly, weekly or daily. We converted these 
primary datasets by making a linear approximation of the growing 
trajectory in order to obtain a mean and a standard deviation by 
month and dataset for each available calendar year. Datasets from 
the same site were assumed to be samples from the same 
population and were not analyzed separately. This could affect the 
value of the whole population mean growth, but here we are only 
interested in the monthly variations of tree growth (site effect is 
considered as a random effect). The final dataset is comprised of 
30 datasets of 25 pantropical forest sites, for a total of 3412 
individual trees (Figure 1, Table 3). The data, monthly tree growth 
mean and SD, are freely available upon request to the 
corresponding author. No specific permits were required for the 
described field studies and this study did not involve endangered 
or protected species. 

Preliminary analysis 

First, we investigated the association between climate variables 
on a monthly time scale through a principal component analysis 
(PCA) on the normalized climate dataset, i.e. climate variables 
were centred and scaled, to describe how the variance of the 
climate dataset was structured. Next, to measure spatial autocor- 
relation in tree growth observations between sites, we computed 
the Moran's Index statistics of the mean monthly growth by sites 
[45]. This index range from -1, strong negative spatial autocor- 
relation, to +1, strong positive spatial autocorrelation. The 
significance of Moran's I is evaluated by using a Z score and p- 
value generated by random permutation [46] . The null hypothesis 
states that there is no spatial autocorrelation for the variable within 
the geographic area. The interannual and intra-annual variability 
climate variable was described by computing the coefficient of 
variation (GV) for each variable (i.e. standard deviation x 100 
divided by the mean) of the annual mean and monthly values of 
the climate time series over the period of the CRU data (1901- 



2009). To detect, estimate and test seasonal patterns in the tree 
growth time series, we used temporal regression models from the 
R package season [47]. The model was fitted using a sine and 
cosine term that together described the sinusoid. These parameters 
were added to a generalized linear model to explain tree growth 
data and test the existence of a seasonal pattern. The existence of a 
seasonal pattern was determined by the zero-test based on 
Snedecor's F statistic. This method is known as the cosinor test. 

Modeling seasonal tree growth 

We modeled the link between tree growth, site and climate 
variables in a mixed linear model framework at a monthly time 
step. Here, we assumed that each growth measurement was 
independent and that our multiple-site design enabled us to deal 
with the collinearity of climate variable present in each site. We set 
the climate variables as fixed effects and the site as a random 
effect. We set the site effect as a random effect in order to avoid 
any statistical bias in our results. This site effect was not further 
analyzed as the main objective of the study was to understand 
seasonal variations of tree growth. In order to rank the climate 
effects, we considered various growth models: 

• mo, the reference model: 

Incij = + C; + e, v (1) 

• wivardim, for each climate variable varclim, an univariate 
model: 

InCij =P 0 + Ci + Pvarclim * Varclmiij + Eij (2) 

• m mc , the model with the best combination of climate 
variables, Combnic, according to the Bayesian Information 
Criterion (BIG): 

Inctj = 

A) + Ci + l } va r clim Xvarclim ij + e -i,i 

varclim^ Comb £ 

Where IriCjj is the average tree diameter growth for the site i for 
the month j, is the random site effect assumed to follow a 
normal distribution J\f(0,o^ jte ) and Bjj is the residual error assumed 
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Figure 2. Principal component analysis of the climate variables. 

dd: cloud cover; pre: precipitation; sol: extraterrestrial solar radiation; 
tmp, tmn and tmx are respectively the daily mean, minimal and maximal 
temperatures; dtr: temperature amplitude; vap: vapour pressure; per: 
potential evapotranspiration; and swc: relative soil water content, (a) 
correlation circle of axis 1 and 2; (b) correlation circle of axis 1 and 3; (c) 
projection of the country classes on the pea axis, India (Indi), Uganda 
(Uga), Brazil (Bra), Mexico (Mex), Costa Rica (CR), Kenya (Ken), Ivory 
Coast (IC), DRC (DRC), Indonesia (Indo), Ethiopia (Eth), French Guiana 



(FG), Ghana (Gha), Ecuador (Ec), Venezuela (Ven), Malaysia (Ma), 
Thailand (Th); (d) projection of the continental classes on the pea axis; 
and (e) projection of the global ecological zones on the pea axis (i.e. 
Tropical rainforest, Tropical moist deciduous forest, Tropical dry forest, 
Tropical mountain system and Subtropical humid forest). Note that axis 
1 and 2 explain 41.70% and 28.12% of the total variation respectively, 
(a). The third axis explained 11.30% of the variance and was linked 
negatively to sol and pre, (b). 
doi:10.1371/journal.pone.0092337.g002 

to follow a normal distribution Af(0,NobSj x a 2 ) weighted by 
1 /Nobsj, which is the inverse of the number of observation periods 
by site. The mixed model parameters were estimated with the 
algorithms of the R package lme4 [48] . To estimate the parameter 
confidence intervals with a probability of 95%, we generated 
10000 parameter samples from the posterior distribution of the 
fitted model parameters using the Markov Chain Monte Carlo 
methods [48]. Then, the confidence intervals for the parameters 
(Highest Posterior Density [HPD] intervals) were constructed from 
the empirical posterior distribution of the 10000 samples as the 
intervals containing the parameters with the nominal probability 
[48]. Evaluation of the model performance was made by 
computing goodness-of-fit-measures with the R package lmmfit 
[49-52]. To find the best variable linear combination that contains 
the maximum of information to link growth and climate variables 
(the model meicj, we ran an exhaustive screening of the candidate 
models using a stepwise procedure based on the Bayesian 
information criterion, BIC [53]. We used BIC, instead of the 
classically used AIC, to avoid over-parameterization and multi- 
collinearity problems, as this criterion is consistent and parsimo- 
nious for model selection with respect to large datasets [54] . We 
made a residual analysis to verify if the error of the model m BIC had 
a bias for any of the selected variables in this model. That is, we 
computed the Pearson's product-moment correlation coeffiecient 
(ft) between the residual of the model m BIC and the variables of this 
model and tested if this coefficient was statistically or different not 
from zero. 

The predictive quality of the fitted models was assessed by 
computing the root mean square errors of predictions, RMSEP. 



RMSEP = 



V™nth ( InCij _ InCiJ f 
Pi pi NstteXNmonth 



(4) 



where IflCu is the observed values of tree growth for the site i 
and for the month j, andj^T. are the model predictions of growth. 

All analyses were performed using the R-project software 
(httpj : / /www. r-proj ec t . org/) . 

Results 

Climate gradients 

The 25 sites represent a large sample of tropical forests under 
different tropical climates corresponding to five global ecological 
tropical zones [55] (Table 4). The gradient of the annual mean of 
precipitation ranges from 973.9 mm (RBSF, Ecuador) per year to 
3948.1 mm (La Selva, Costa Rica). The coefficient of variation 
(CV) of interannual precipitation ranges from 12.5% (SERS, 
Thailand) to 35.2% (RBSF, Ecuador). Intra-annual CV of 
monthly precipitation distribution precipitation ranges from 13% 
(CPM,Brazil) to 116.9% (El Palmar and La Barcinera, Mexico). 
The sites undergo a large variability in the distribution of 
precipitations throughout the year, from zero months (La Selva) 
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Table 5. Seasonality of mean monthly tree growth with cosinor analysis for the studied sites. 



Site 


Amplitude 


Phase 






Low phase 


p value 


Attapadi 


0.129 


Month 


= 


7.6 


Month 




1.6 


<0.05 


Budongo 


0.031 


Month 


= 


7.3 


Month 




1.3 


<0.05 


CPM 


0.104 


Month 


= 


12.2 


Month 




6.2 


<0.05 


El Palmar 


0.437 


Month 


= 


7.7 


Month 




1.7 


<0.05 


FLONA SFP 


0.213 


Month 


= 


1.1 


Month 




7.1 


<0.05 


Guanacaste 


0.396 


Month 


= 


7.3 


Month 




1.3 


<0.05 


Ibicatu 


0.470 


Month 


= 


1.5 


Month 




7.5 


<0.05 


Kakamega 


0.019 


Month 


= 


7.5 


Month 




1.5 


<0.05 


La Barcinera 


0.361 


Month 


= 


6.8 


Month 




12.8 


<0.05 


La Selva 


0.092 


Month 


= 


9 


Month 


= 


3 


<0.05 


Lamto 


0.067 


Month 


= 


6.9 


Month 




12.9 


<0.05 


Luki forest 


0.535 


Month 


= 


1.4 


Month 


= 


7.4 


<0.05 


Marajoara 


0.684 


Month 


= 


12.6 


Month 




6.6 


<0.05 


Muara Bungo 


0.211 


Month 


= 


8.6 


Month 


= 


2.6 


<0.05 


Munessa-Shashamene Forest 


0.380 


Month 


= 


6.4 


Month 




12.4 


<0.05 


Paracou 


0.091 


Month 


= 


2.6 


Month 




8.6 


<0.05 


Pinkwae 


1.033 


Month 


= 


7 


Month 




13 


<0.05 


RBSF 


0.361 


Month 


= 


2.2 


Month 




8.2 


<0.05 


RFC 


0.645 


Month 




7.6 


Month 




1.6 


<0.05 


Rio Cachoiera 


1.235 


Month 




1.3 


Month 




7.3 


<0.05 


Selangor plantation 


0.380 


Month 




6.7 


Month 




12.7 


<0.05 


SERS 


0.033 


Month 




7.4 


Month 




1.4 


<0.05 


Tapajos 


0.089 


Month 




1.9 


Month 




7.9 


<0.05 


Tinte Bepo 


0.223 


Month 




6.9 


Month 




12.9 


<0.05 


ZF-2 


0.066 


Month 




1.1 


Month 




7.1 


<0.05 



Month = 1 corresponds to January. The amplitude is the difference between the higher and lower points in mm.month -1 of the sinusoid fitted in the cosinor analysis. 
The phase and the low phase are, respectively, the month with the highest/lowest tree growth value according to the sinusoid fitted in the cosinor analysis. A 
Pvalue<0.05 indicates that the statistically significant existence of a seasonal pattern cannot be rejected. 
doi:1 0.1 371 /joumal.pone.0092337.t005 



to 8-9 months with less than 100 mm of rain (Lamto, Tinte Bepo, 
Pinkwae, RBSF). The relative soil water content exhibits the same 
intra-annual pattern, ranging from sites without seasonality, intra- 
annual CV SWC < 3.5% (Cpm, La Selva, Muara Bungo, Rio 
Cachoiera and Selangor), to sites with an intra-annual CV„ > 
37.8% (El Palmar, la Barcinera, RBSF). The annual average of 
relative soil water content is above 60% for 24 of the 25 sites. The 
annual cloud cover mean ranges from 42.9% (Munessa-Shasha- 
mene, Ethiopia) to 87.2% (Luki Forest, DRC) and its intra-annual 
variation is similar to the intra-annual variation of precipitation. 

Three major climate gradients are represented in the Principal 
component analysis (PCA)(Figure 2). The first axis (x-axis in 
Figure 2a) represents a gradient of increasing temperature, that is 
highly correlated with vapour pressure. Sites close to the equator 
that are located at high altitudes are on the left of this axis, with 
sites far from the equator experiencing a cold season (Figure 2e). 
The second axis (y-axis in Figure 2a) has a strong contribution 
from variables that reflect water availability. The third axis (y-axis 
in Figure 2b) appears to be related to solar radiation. Ordination 
of the continental location of the studied sites on the PCA axis 
(Figure 2d) revealed that sites in Asia were slightly warmer and had 
higher vapour pressures than those in America and Africa. 
Ordination of the studied countries on the PCA axis (Figure 2c) 
reveals that some countries are separated from others because of 



climate: from cold (Ecuador) to warm (Indonesia) and from dry 
tropical (Mexico) to wet tropical (Malaysia)(Figure 2e). In the 
following analysis, we kept all the climate variables as we have a 
strong physiological assumption of their effect on tree growth 
(Table 1). 

Tree growth descriptive analysis 

We cannot reject the null hypothesis that there is no spatial 
autocorrelation present in mean tree growth observations between 
sites at alpha = 0.05 (Moran's I observed = 0.0102, expected = — 
0.042, SD 0.123, p.value = 0.672). Tree diameter growth shows a 
significant intra-annual seasonality at all the studied sites, even 
sites with mean annual precipitation rates close to 4000 mm (La 
Selva, Costa Rica)(Table 5 and Figure 3). 

Climate effects on tree growth 

More than 28.69% of the observed seasonal variation in tree 
growth may be imputable to the site effect (reference model m 0 , 
Table 6), while climate variables alone explain a maximum of 
19.82% (squared root monthly precipitation). The variables sol, eld, 
vap, dtr, tmn and swe explained between 9.65 and 16.30% of the 
climate effect, while tmp, pet and tmx explained less than 2.13%. 
The selection procedure, which used the BIC criterion, kept four 
climatic variables [Jpre, sol, dtr, swe) in the final multivariate 
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Figure 3. Evolution of mean monthly tree growth values and standard deviation for the studied forest sites. The red line represents a 

cubic smoothing spline. 

doi:10.1371/journal.pone.0092337.g003 



model m Blc (Table 6). These four climate variables together 
explained 29.79% of the total observed monthly tree growth 
variation. The parameter values for the fixed effect in the 
univariate analysis indicate the direction of the relation between 
the climate variables and tree growth (Table 6). Among the four 
selected variables, pre (0.03), swc (0.58) and sol (0.09) have a positive 
link with tree growth and dtr (—0.08) has a negative link (Table 6). 
The obtained root mean squared error of prediction (RMSEP) 
was slighdy below the mean value of observed monthly 



growth in diameter (mean growth 



0.325 mm.monfh 



mated growth when it was above 1 mm.month - (Figure 4). 
There was no significant correlation between the selected variable 
and the residuals of the model m BIC (p = 4.6 x 10 -6 , pvalue > 
0.99; p sal = 1.8 xl0~ 5 , pvalue > 0.99; p itr = 1.9 xlO" 4 , pvalue > 
0.99; p m , c = 6.1x10 6 , pvalue > 0.99). Tree growth values 
increase linearly with extrasolar radiations (Figure 5), and a strong 
increase in tree growth occurred between 0 and 200 mm.month of 
precipitation; above these values, the increase in tree growth was 
less pronounced and had a linear shape. 



RMSEP = 0.279 mm.month ). In general, the model underesti- 
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Figure 4. Observed versus predicted diameter at breast height 
(dbh) growth under the model m Blc . The red line is the identity line 
y = x. Note that the model underestimated the diameter growth above 
1 mm.monthT 1 . 

doi:1 0.1 371 /journal.pone.0092337.g004 




'On, 



200 1<c ,v 



Figure 5. Predicted diameter at breast height (dbh) growth 
under the model m BIC and relationship with precipitation and 
extra solar radiation. Predictions were computed using pre, sol, 
mean dtr, mean swc and mean of the random effect. The extra solar 
radiation unit is equivalent of evaporation in mm. month -1 , 
1 mm.month -1 equivalent of evaporation = 2.45 MJ.m~ 2 .month _1 . 
doi:10.1371/journal.pone.0092337.g005 



Discussion 

Our meta-analysis provides the first overview of climate effects 
on tree growth variations at a seasonal time scale for tropical forest 
areas. Tree growth reveals a strong intra-annual seasonality at all 
sites (Table 5), including seasonally dry sites and wet tropical 
forests. We found a common climatic signal explaining 29.79% of 
the observed seasonality in forest growth (Table 6). Precipitation 
(19.82%) and extraterrestrial solar radiation (16.30%) were the 
major climate drivers. The tree growth average per site (i.e. the 
random site effect) explained 28.69% of overall growth variation, 
and a substantial fraction of variation (50%) in growth remained 
unexplained with our modeling framework. We have to acknowl- 
edge than our sample of 3412 trees could not be used to fully 
reproduce the complete behavior of all the tropical forests, 
however our analysis demonstrated that a common signal in the 
climate determinism of tree growth seasonality is observed at the 
pan-tropical scale. 

Climate effect on tree growth 

In this study, precipitation was found to be the major seasonal 
driver of tree growth. Precipitation strongly impacts tree growth, 
as directly observed in seasonal and unseasonal tropical climates 
[19,35,56-59] and as deduced from experimental forest droughts 
[60,61]. The prominence of precipitation as a predictor of forest 
growth is slightly different from previous studies [62]; in French 
Guiana, we found that soil water availability was the main 
determinant of diameter growth, i.e. better than rainfall. Of 
course, both variables are strongly correlated (Figure 2), but this 
difference could be explained by some sites where relative soil 
water availability swc was at its maximum throughout the year 
despite marked precipitation seasonality (e.g. La Selva, CPM, 
Muara Bungo, Rio Cachoiera and Selangor). At the same time, 
the importance of solar radiation, sol, in the complete model (niBjcj 
reflects the obvious role of light in shaping tree growth (Figure 5). 



Solar radiation is directly linked to PPFD (Photosynthetic Photon 
Flux Density), which in turn drives carbon uptake and plant 
growth [29] . Some authors support the hypothesis that increasing 
surface solar radiation contributes to the increasing forest growth 
rate over the Amazon [63,64]. It must be noted that the 
extraterrestrial solar radiation we used has a value above the real 
solar radiation value reaching the forest surface. Indeed, sol was 
computed as a monthly mean over 50 years and does not account 
for local cloud cover or aerosol radiation absorption. 

Investigating the effects of temperature on the physiology of 
tropical forest trees [31,65] is of primary importance today given 
the temperature increases expected over the next century [5,7]. 
Some authors suggest that tropical trees are more sensitive to 
temperatures than other trees because (i) they live at or close to the 
highest annual average temperatures on Earth, and (ii) tropical 
species naturally encounter limited variation in temperatures (< 
4°C over 20° of latitude) [32]. Our results suggest that temperature 
variations are of secondary importance in shaping tropical tree 
seasonal growth; nevertheless, they do play a role. Minimal 
temperature was slightly positively correlated with tree growth, 
whereas maximal temperature had no effect. This positive relation 
between tree growth and temperature is not consistent with 
previous observations in Costa Rica [35,66], where an increase in 
night-time temperatures had a negative effect on tree growth, and 
in East Africa [19], where maximum temperatures had a negative 
effect on tree growth. Such patterns were not found at our global 
working scale. 

Tree growth variability in time and space 

Our results suggest that most tropical trees experience seasonal 
cycling growth even in extremely wet environments, as already 
reported at la Selva [35]. On a pantropical scale, we showed that 
secondary growth was higher during the wet season; tree growth 
increased with precipitation and relative soil water content 
(Table 6). This result is not consistent with other studies that 
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have shown that tropical forests are able to maintain or even to 
increase their productivity during the dry season [67-70]. This 
suggests different uses of carbohydrates. Some studies have 
reported that evergreen species in seasonally dry environments 
accumulate carbohydrates during the dry season because photo- 
synthesis continues while wood production ceases [71] and that 
deciduous species accumulate carbohydrates at the onset of the dry 
season to support respiration costs when they are leafless [72,73]. 
Some studies relied on the Enhanced Vegetation Index (EVI, an 
index of canopy photosynthetic capacity [30]) to highlight changes 
in forest phenology driven by the solar cycle [74,75]. In the same 
way, a variation in the EVI has been observed for the entire 
Amazon region [30], and the link between phenology and the 
solar cycle has also been reported for Terra Firme forests [76]. 
The variability in tree growth at a seasonal scale is likely driven by 
climate seasonality and dependent on the seasonal allocation of 
carbohydrates to processes other than tree diameter growth, such 
as leaf and root production or respiration. 

Stem shrinkage during dry periods may be an important 
limitation of this work [77-79], e.g. the negative value of mean 
annual tree growth at Pinkwae. Other negative monthly growth 
values exist at almost all the study sites. In a tropical forest in 
Ethiopia experiencing a strong seasonality, high-resolution elec- 
tronic dendrometers have been combined with wood anatomy 
investigation to describe cambial growth dynamics [80]. These 
authors concluded that water scarcity during the long dry season 
induced cambial dormancy. Furthermore, after the onset of the 
rainy season, (i) bark swelling started quite synchronously among 
trees; (ii) bark swelling was maximum after few rainy days; and (iii) 
evergreen trees were able to quickly initiate wood formation. 
Recently at the Paracou forest site, some authors have showed that 
biomass increments were highly correlated between the first and 
the last quantile of trunk bark thickness and between the first and 
the last quantile of trunk bark density, suggesting that secondary 
growth is driven by cambial activity [75]. 

In this study, we focused on seasonal variation of tree growth, 
but the inter-site variance must be studied to build a full predictive 
model. At our pantropical working scale, there was no evident 
spatial auto-correlation. The site effect, i.e. the average growth of 
trees at a given site, is likely shaped by several additional 
environmental variables, such as soil fertility, forest floristic 
composition and competition for light and nutrients [81-83]. 
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Tropical tree growth under climate change 

Globally, current IPCC scenarios predict an intensification of 
the dry season in tropical areas during the 21st century [7]. 
Amongst climate variables, our results highlighted the predomi- 
nant role of precipitation in shaping seasonal forest growth. If a 
global rainfall reduction is confirmed in the future, it can be 
expected that tree growth will be affected (Figure 5). Climate 
change is not the only possible cause for reductions in 
precipitation; deforestation could lead to reductions in precipita- 
tion as the air passages over forests increase tropical rainfall 
intensity [4]. Solar radiation appeared almost as important as 
rainfall. Cloud reduction due to drier dry seasons, may 
subsequently improve tree growth. Our results suggest that, with 
the global changes observed, forest productivity due to tropical 
tree growth will encounter modifications due to change in climate 
seasonality. 
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